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Abstract 

We have calculated differential, total and transport cross sections 
for muonic hydrogen in the n = 2 state scattering from hydrogen. 
The metastable fraction of the 2S state that slows down below the 
2P threshold without undergoing collisional quenching has been cal- 
culated as a function of the initial kinetic energy using a Monte Carlo 
kinetics program. Contrary to earlier estimates, the metastable frac- 
tion in the kinetic energy range of 2 — 5 eV cannot be neglected. 



1 Introduction 



The 25 state of muonic hydrogen offers interesting possibilities to do precision 
tests of QED and to determine the proton RMS charge radius (see |l| and 
references therein). An isolated {fJLp)%s is metastable with a lifetime mainly 
determined by muon decay (about 2.2 /xs). In liquid or gaseous hydrogen the 
lifetime of the 25 state is shortened considerably because of Stark mixing 
followed by 2P —> IS radiative transitions. If a sizeable fraction of muonic 
hydrogen atoms ends up in the 2S state with a sufficiently long lifetime, then 
precision laser experiments with this metastable 2S state become feasible. If 
the (/xp)25 has kinetic energy below the 2P threshold (laboratory kinetic 
energy T = 0.3 eV), then Stark transitions 25 — ► 2P are energetically 
forbiddenP]. The metastable fraction of (np)2S m hydrogen depends on the 
kinetic energy at the time of formation. 

The first estimate of the (//p)2S lifetime was done by Kodosky and Leon 
H. They calculated the inelastic 25' — > 2P cross section in a semi classical 
framework and concluded that the 25 state for T > T will be rapidly depop- 
ulated except for very small target densities. However, this model did not 
consider deceleration due to elastic 25 — > 25 scattering. A more elaborate 
approach was developed by Carboni and Fiorentini || . They calculated both 
elastic 25 — ► 25 and inelastic 25 — ► 2P cross sections quantum mechanically 
and estimated the probability for a (fip)2S atom to slow down below thresh- 
old from a given initial energy. The results of their calculations show that 
a sizeable fraction of (fip)2s formed at kinetic energies less than 1.3 eV can 
slow down below the 2P threshold. 

The metastable fraction of (//p)2S per stopped muon can in principle 
be calculated in a cascade model which takes the different processes (Stark 
mixing, radiative decays, etc.) into account [[7],|8|. However, if one knows the 
fraction of stopped muons which reaches the 25 state (regardless of energy) 
and the kinetic energy distribution on arrival in this state, then it is sufficient 
to treat the final part of the cascade (n = 1,2). This information can be 
obtained from experiments. The fraction of stopped muons which arrives 
in the 25 state can be determined from the radiative yields @-ll|: it was 



found in Ref. that between 2% and 7% of the fip reach the 25 state 
in the pressure range 0.33 — 800 hPa. The kinetic energy distribution for 
HP in the 15 state, which for low pressures is expected to be very similar 
to that of the 25 state just after arrival, can be obtained from diffusion 



experiments [12, 13|. The median energy is found to be about 1.5 eV for a 



Some quenching will, however, occur because 2S — 2P mixing during collisions allows 
radiative transitions to the IS state (See Refs. [0-0). 
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Figure 1: Coordinates used in the calculations: R is the vector from the 
target proton to the center of mass of the up, r is the relative vector of the 
fip system. 



target pressure of 0.25 hPa |13 | 



The purpose of this paper is to calculate the fraction of fip in the 25 state 
which reaches kinetic energies below the 2P threshold as a function of the 
initial kinetic energy T. We will also present a fully quantum mechanical 
calculation of fip+H. differential cross sections which are used in our Monte 
Carlo simulation of the kinetics. 

The paper is organized as follows. The theoretical framework of the 
quantum mechanical calculation of the cross sections is outlined in Section 2. 
The calculated cross sections are discussed in Section 3. Section 4 presents 
the calculations of the metastable 25 fraction. The summary of the results 
is given in Section 5. 

Unless otherwise stated, atomic units (h = a>o = m e = 1) are used 
throughout this paper. The unit of cross section is Oq = 2.8 • 10 _17 cm 2 . 



2 Quantum Mechanical Approach to the Cal- 
culation of the Cross Sections (fip) n i + H —> 

{W)nl> + H 

For the benefit of the reader, we briefly describe the quantum mechanical 
calculation of fip+B. scattering in the coupled-channel approximation. The 
three-body wave function ip{r, R) where the coordinates r and R are defined 
in Fig. |l| satisfies the Schrodinger equation 

Ht/j(r,R) = Et/j(r,R) (1) 
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where the Hamiltonian is given by 



H = ~+H„ + V{r,R) . (2) 

Here fi = m p m^ p / (m p + m MP ) is the reduced mass of the p — ftp system, 
with m p being the proton mass and ra w the total fip mass. The two-body 
Hamiltonian of the fip atom, H W) includes the Coulomb interaction and a 
term that describes the shift of the nS state (mainly because of the vacuum 
polarization) with respect to the states with I > 0. For the case n = 2 
considered below, the 25* state is lower than the IP by AE = 0.21 eV. 
The much smaller fine and hyperfine structure splitting is neglected. The 
potential V(r, R) describes the interaction of the \ip system with the target 
protonQ: 

^-■^Ijfbr ia+ti-eH (3) 

where e = m^/m^p = 0.101. 

Equation ([!]) is solved in the coupled-channel approximation by using a 
finite number of basis functions to describe the state of the fip. For the 
problem of nlm — > nl'm! scattering considered in this paper, the set of n 2 
eigenstates with principal quantum number n has been selected but the basis 
can be extended in a straightforward manner. With n fixed, let Xim{ r ) denote 
the normalized eigenf unctions of the atomic Hamiltonian H w with the energy 
E n i, the square of the fip internal angular momentum 1 (eigenvalue 1(1 + 1)) 
and its projection along the z-axis l z (eigenvalue m). The total wave function 
ip(r,R) is expanded as follows 

iKr, R) = R- 1 £ tMRWin, r) (4) 

JMLl 

where 

y J L f{n,r) = (LlM L m\JM)Y LML (fl) X i m (r), ^ = R/R . (5) 

are simultaneous eigenfunctions of J 2 , L 2 , l 2 and J z with eigenvalues J(J+1), 
L(L + 1), 1(1 + 1) and M, respectively. Here L is the p — fip relative angular 
momentum, J = L + 1 is the total orbital angular momentum of the system. 

2 For the sake of simplicity, we ignore the fact that the protons are identical particles. 
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For a given value of J the system of radial Schrodinger equations has the 
form 

+ = o 

(6) 

where the potential matrix elements are calculated in the basis @: 

(n,r\LUM)=y J J I (n,r) . (7) 

The matrix elements of the potential (|3p have been calculated analytically; 
the corresponding formulas are rather lengthy and will be given elsewhere. 
From the asymptotic form of the solution of the n 2 coupledf] equations (^|), 
the scattering matrix S is extracted and cross sections can be calculated 
using standard formulas. The scattering amplitude for nlm — > nl'm' is given 
by 

4:71 

fnim->ni>m>(n) = — j== £ i L ~ L Y l <m' l («) {L'l'M' L m'\S - 1 \LlOm)Yl (0) 

(8) 

where 

(L'l'M' L m'\S\LlOm) = Y,(L'l'M' L m'\JM)(JM\LlOm)(L'l'J\S\LU) . (9) 

JM 

As a consequence of rotational symmetry, the matrix elements (L'l'J\S\LU) 
do not depend on the quantum number M. The differential cross sections 
for the transitions nl — > nl' are given by 

dn ~ (2i + i) v k Unl " l ' m '\ 1 j 

V / mm 

where k and k' are the magnitudes of the relative momenta in the initial and 
final state, correspondingly. 

The total cross sections of the transitions nl — > nl' have the form 



3 Because of parity conservation the equations decouple into two sets of respectively 
n(n + l)/2 and n(n — l)/2 coupled equations. 
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and the corresponding transport cross sections are given by 



^ = /^(l-cos0)^|^ . (12) 

In order to treat the long distance behaviour of the /ip+H interaction 
properly, the effect of electron screening must be taken into account. This 
is done by multiplying the matrix elements (L'l'JM\V\LlJM) in Eq. by 
the screening factor 

F(R) = (l + 2R + 2R 2 )e~ 2R (13) 

which corresponds to the assumption that the electron of the hydrogen atom 
remains unaffected in the IS state during the collision. 

For p — [ip separations R smaller than a few units of the up Bohr ra- 
dius, a M = 0.0054, our model cannot be expected to be valid because the 
truncated set of basis functions in Eq. ([|) is not sufficient to describe the 
total three-body wave function ip(r,R). Furthermore, exchange symmetry 
between the two protons must be taken into account. We can estimate the 
sensitivity of our results to the short range part of the interaction by using 
the dipole approximation for the potential (|3|). The interaction in the dipole 
approximation is given by the first nonzero term in the expansion of Eq. (|3]) 
in inverse powers of R: 

Vda{t,R) = ^ = ^ • (14) 

A certain problem arises in the dipole approximation for a few low partial 
waves ( J < 5): the Schrodinger equation becomes ill defined because of the 
attractive 1/R 2 singularity^. Following Ref. we cure this difficulty by 
placing an infinitely repulsive sphere of radius r m j n around the target proton. 
The sensitivity of the results to this cutoff parameter r min will be used below 
as an estimate of the importance of detailed description of the interaction at 
short distances. 

In this paper we are interested in the 21 — > 21' transitions, and only four 
states n = 2 are used to describe the \ip part of the total wave function. 
The four coupled second order equations @ are solved numerically for J = 
0, 1, J max where the highest partial wave J max is chosen large enough to 
ensure the convergence of the partial wave expansion at given collision energy. 



4 This is a problem only in the dipole approximation. The exact matrix elements are 
all finite for R = 0. 
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Until now we have considered the HP collisions with the atomic target. 
Treating the collisions with hydrogen molecules is a formidable task (even 
for fip in the ground state [Q) which we do not attempt here. The inelas- 
tic threshold To for 2S — > 2P transitions is 0.44 eV for atomic target and 
0.33 eV for a molecular target. To get the correct threshold value for the 
inelastic cross sections one can substitute the atomic hydrogen mass with 
the molecular one. By varying r m ; n and the target mass one can obtain an 
estimate of the theoretical uncertainty of our approach. 

The present model for calculating cross sections for (/zp) n=2 — H scattering 
is a straightforward extension of the one by Carboni and Fiorentini . There 
are three major differences: we solve the four coupled differential equations 
exactly while Ref. |J treated non adiabatic terms as a perturbation. The 
second difference is that we include the full angular coupling while Ref. 
omitted some minor terms. Finally, we use exact matrix elements for the 
HP — p interaction while Ref. considers only the dipole approximation. 
These approximations were justified by Ref. M as follows: for small kinetic 
energies the velocity of the muonic hydrogen atom is so low that the motion 
can be regarded as nearly adiabatic. The angular coupling terms that were 
omitted in Ref. are of the order 1 which is much smaller than the remaining 
angular coupling terms of the order of J(J+ 1) (the angular momenta as high 
as J ~ 15 contribute to the 25* — > 2P cross section at T = 1 eV). The electric 
field from a hydrogen atom is strong enough to induce Stark transitions in 
the up for the distances R ~ a . Therefore, the regions where the dipole 
approximation is valid (R 3> a^) are supposed to be most important. The 
diffusion experiments |12|,p^] have shown that a sizeable fraction of (fip)n=2 
atoms has kinetic energies of several eV. In this high energy region the non 
adiabatic couplings become strong and the model of Ref. can not be 
expected to give accurate results. 

The problem of /ip+H scattering has been treated fully quantum me- 
chanically in Refs. |I6|-[T9|1. However, these calculations did not include the 
25 — 2P energy splitting and the question of the metastability of (hp)2S was 
not addressed. Stark mixing has been studied in the semiclassical straight- 



line-trajectory approximation in Refs. I5[ [2(| . We have calculated the 
2S — > 2P Stark mixing cross sections in the semiclassical approach in order 
to compare with the quantum mechanical results. A more detailed compar- 
ison between semiclassical and quantum mechanical calculations of HP + H 
scattering will be given elsewhere. 
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Figure 2: Transport cross section 25 — > 25 (solid lines) and Stark mixing 
cross section 25* — > 2P (long-dashed lines) vs. laboratory kinetic energy T. 
The thick lines are the results of the present calculations, the thin lines are 
obtained from Fig. 3 of Ref. || and the dotted line shows the result of the 
semiclassical calculation in the straight-line-trajectory approximation. The 
quantum mechanical cross sections are calculated with M targ et = Mh 2 . 

3 The Cross Sections of (fip)2i Scattering from 
Hydrogen 

Using the method described in Section 2 the S-matrix has been calculated 
for the laboratory kinetic energy range Tq < T < 6 eV. Unless otherwise 
explicitly stated, the results shown are obtained with the exact potential (^|). 
Electron screening is always taken into account. Both atomic and molecular 
mass of the target (Mt arget = Mr, Mh 2 ) have been used. 

Figure ^| shows the 25 — > 25 transport cross section and the 2S — > 2P 
Stark mixing cross section in comparison with the results from Ref. || (the 
molecular mass is used in both cases). There is a good agreement for the 
Stark mixing 25 — > 2P cross section below 1.7 eV. For the 25 — > 25 transport 
cross section, the agreement is fair, with the discrepancy being typically less 
than 30%. 

In order to estimate the theoretical uncertainty of our approach, we cal- 
culated the cross sections in the dipole approximation with short distance 
cutoff 0.01 < r min < 0.05 for M target = M H and M target = M Ha . For fixed 
value of M target , the cross sections for the three reactions 25 — > 2P, 2P — > 25 
and 2P — > 2P are weakly dependent on r m j n . This shows that these reactions 
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Figure 3: Partial wave cross sections for the reactions 25' — > 2P (a) and 
2S — > 2S (b) at laboratory kinetic energy 3 eV and M target = M H . The 
histograms show the results for exact matrix elements; the dipole approxi- 
mation with r min = 0.01 and r min = 0.05 is shown by o and x, respectively. 
The semiclassical result for the 2S — > 2P transition is shown with a solid 
line. 



are dominated by the long range part of the interaction V(r, R). The only 
process rather sensitive to the value of r min is the elastic scattering 2S — > 2S. 
This can be understood by considering the adiabatic energy curves for low 
angular momentum. The energy curve which corresponds asymptotically to 
the 2S state is attractive while those corresponding to the three 2P states 
are repulsive. Therefore, in the adiabatic approximation the 25 — > 2S cross 
sections are expected to depend on the short range part of the potential 
while this is not the case for 2P — > 2P scattering. At energies above 2 eV 
the semiclassical approximation is in a good agreement with our quantum 
mechanical results. However, this semiclassical approximation does not treat 
the threshold behaviour correctly. 

A more detailed comparison of the approximations used can be done by 
plotting the J dependence of the partial wave cross sections <jj at fixed energy 
as shown in Fig. |3|. The quantum mechanical results obtained for the exact 
potential and the dipole approximation with the short range cutoff agree well 
for angular momentum J > 5 while the lowest partial waves are sensitive to 
the short range behaviour of the approximating potentials. The reason is 
that for J > 5 the centrifugal barrier is strong enough to prevent the fip 
from approaching close to the target proton. For J < 5 the up can get very 
close to the proton and the use of a small number of atomic orbitals is not 
sufficient — a better description in this region is needed in a true three- 
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Figure 4: Differential 25* — > 2S cross sections vs. CMS scattering angle 9 for 
three different laboratory kinetic energies: (a) M targct = M H , (b) M targct = 
M H2 . 



body framework. It is seen that a substantial part of the 2S — > 2S cross 
section comes from partial waves with low J, so this result also explains why 
the uncertainty of the elastic 2S cross section is larger than for the other 
reactions. The semiclassical calculation can be compared with the partial 
wave cross sections by using the relation between the impact parameter p, 
the relative momentum k and the angular momentum J 

kp = J + 1/2 . (15) 

For large J (large impact parameter) there is a very good agreement between 
the semiclassical contribution to the 2S — > 2P cross section and the quantum 
mechanical partial wave result. An example of the differential cross sections 
for the reaction 2S — > 2S given in Fig. 4 shows a characteristic pattern with a 
strong forward peak and a set of maxima and minima, which is in qualitative 
agreement with Ref. [[Tj| where the adiabatic approach was used. 



4 The Surviving Fraction of the Metastable 
(Mp)25 State 

The surviving metastable fraction f(T) is defined as the probability that 
the up atom in the 25 state with initial kinetic energy T reaches the energy 
below the 2P threshold by slowing down in elastic collisions. Assuming that 
the rate of the radiative transition 2P —* IS, \2P^is = 1-2- 10 11 s _1 , is much 
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Figure 5: Collisional rates for 21 — > 2Z' in liquid hydrogen (N — N , M target = 
Mr) and the radiative 2P — > IS 1 transition rate. 



larger than the Stark mixing ratef], the surviving fraction /(T) was estimated 
in || by the formula 

= exp M™/* + ;fW) 2 ^ a 2 ^2P(T0 \ g) 



with M target = M H2 . It was found that a sizeable fraction of (fip)2s atoms 
formed at kinetic energies below 1.3 eV slows down below threshold. 

Equation fllEp is based on the approximation of continuous energy loss. 
To provide a more realistic treatment of the evolution in kinetic energy we 
use a Monte Carlo program based on the differential cross sections for the 
four processes 2S -> 2S, 2S -> 2P, 2P -> 25 and 2P -> 2P. In addition to 
the collisional processes, the 2P — > IS 1 radiative transition is also included 
in the code. The fate of a fip formed in the 2S state with kinetic energy T is 
thus either to undergo 2P — > IS radiative transition after the Stark mixing 
25 — > 2P or to end up in the 2S state with kinetic energy below the threshold 
with probability f(T). Figure | shows the rates, \2i-,2i' = N v<j2i-*2i' , for the 
collisional transitions in liquid hydrogen in comparison with the radiative 
deexcitation rate \2p^is- In liquid hydrogen the Stark mixing rates are so 
large that the /ip states are expected to be statistically populated for kinetic 
energies T > 2 eV (where threshold effects can be neglected). 

5 With our result for the Stark mixing rate 2P — > 25 at 1 eV as a function of the target 
density N, \2P^2S = Nv<i2P^2S ~ 4 • 10 12 (iV/7Vo) s _1 where Nq is the liquid hydrogen 
density 4.25 • 10 22 atoms/cm 3 , the range of validity is N <§C 0.03iVo. 
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Figure 6: The metastable surviving fraction f(T) of the (fip)2S states vs. 
initial kinetic energy T. The thick and thin solid lines show the results of 
Monte Carlo calculations with M target = Mh and M tar get = Mh 2 , respectively. 
The dotted line shows the result of Carboni and Fiorentini and the long- 
dashed line shows f(T) calculated using the method of Ref. || (Eq. (|16|)) 
with our cross sections. The results are valid for densities iV < 10~ 2 iVo. 



Figure |] shows the surviving fraction f(T) calculated with the Monte 
Carlo program for target density KT 6 < N/N < KT 2 . The approxima- 
tion (lC) gives somewhat higher values for the survival probability than the 
exact kinetics calculation at T < 1.4 eV. The Monte Carlo results at high en- 
ergies (T > 1.5 eV) are significantly larger than those obtained from Eq. ( |16| ) 
where continuous energy loss is assumed. The reason is that the backward 
scattering (see Fig. 4) with maximum possible energy loss plays an important 
role in bringing the {np)2S atoms below the 2P threshold for higher energies. 

In order to estimate the theoretical uncertainty of f(T) we performed the 
Monte Carlo calculation with the cross sections obtained in the dipole ap- 
proximation. In all cases the Monte Carlo calculations are consistent with the 
results corresponding to the exact potential: at 2 eV the surviving fraction is 
in the range 15 — 20% for atomic target mass and 10 — 16% for molecular tar- 
get mass. The use of the target mass Mh instead of Mh 2 leads to somewhat 
higher survival fractions because of a simple kinematical reason: the loss of 
kinetic energy in a collision with the same angle in the CMS is larger for 
the target of smaller mass and, furthermore, the inelastic threshold is higher 
for the scattering from the atomic target. Merely substituting the atomic 
mass with the molecular mass for the hydrogen target does not account for 
the additional energy loss due to rotational and vibrational excitations of 
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H 2 . One would therefore expect that the slowing down process is more effi- 
cient than this model suggests and the survival probability calculated with 
molecular target mass is underestimated. The opposite is true for calcula- 
tions with atomic target: here the transfer of kinetic energy from the fip to 
the individual hydrogen atoms is not restricted by molecular bindings. Thus 
results with atomic target probably give somewhat optimistic results for the 
surviving fraction. 



5 Conclusion 

The main results of this paper can be summarized as follows. The detailed 
Monte Carlo kinetics calculations predict the surviving metastable fraction 
of the 2S state of up to be larger than 50% for the initial kinetic energy 1 eV 
in agreement with earlier estimates ||. For higher initial kinetic energies, 
our result is significantly larger than the earlier estimates: the surviving 
metastable fraction for T = 5 eV is about 4%. This effect is due to a sizeable 
contribution of backward scattering in elastic collisions. 

Our Monte Carlo calculations are based on the cross sections calculated in 
the coupled-channel approximation. The main limitation of this method for 
the problem concerned comes from the use of a small number of atomic states 
to describe the \ip system and the neglect of the molecular structure of the 
target. A more accurate treatment of the fipp three body problem is needed 
in order to do reliable calculations for a few lowest partial wave amplitudes. 
Our approach, however, is well suited for the description of the collisions 
with the characteristic scale of impact parameters of the order of which 
is exactly the case for the problem involved. Therefore our results provide 
a significantly improved basis for a better estimate of the metastable {np)2s 
fraction [2TJ which is very important for the planned Lamb-shift experiment 
at PSI 

Further details and more results concerning the scattering of the [ip atoms 
in the excited states n > 2 will be published elsewhere. 
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